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Antiferromagnetism of ultracold fermions in an optical lattice can be detected by Bragg diffraction 
of light, in analogy to the diffraction of neutrons from solid state materials. A finite sublattice 
magnetization will lead to a Bragg peak from the {\\\) crystal plane with an intensity depending 
on details of the atomic states, the frequency and polarization of the probe beam, the direction and 
magnitude of the sublattice magnetization, and the finite optical density of the sample. Accounting 
for these effects we make quantitative predictions about the scattering intensity and find that with 
experimentally feasible parameters the signal can be readily measured with a CCD camera or a 
photodiode and used to detect antiferromagnetic order. 

PACS numbers: 37.10.Jk, 42.50.Ct, 71.10.Fd 



I. INTRODUCTION 

Ultracold atomic gases in optical lattices may be used 
as analog quantum simulators of condensed matter mod- 
els Simulating the Fermi-Hubbard model is particu- 
larly interesting. This model is nearly impossible to solve 
using traditional techniques, yet is amenable to simu- 
lation with cold atoms. Most importantly, this model 
is relevant for understanding a wide class of strongly 
correlated electron systems: for example, many believe 
that it captures the essential physics responsible for high- 
temperature superconductivity in cuprates [2J. The ma- 
jor technical challenges to this program are developing 
improved cooling and measurement techniques. Here we 
explore how light scattering can be used to detect "mag- 
netic" order in a gas of fermionic atoms in an optical 
lattices, focussing on the antiferromagnetic order found 
at low temperatures in the insulating phase of the Hub- 
bard model. Producing (and detecting) such an antifer- 
romagnetic state is a key step on the road to exploring 
superfluidity in the Hubbard model. 

The Hubbard Hamiltonian for a two-spin-componcnt 
Fermi gas is 

H = -t ^ (4,* c j,<? + h - c -) + u X] ni >t n U> 

<»,j>.o-=T,i « 

where i,j label lattice sites, angle brackets denote sums 
over nearest neighbors, a is the spin label, and the oc- 
cupation of site i is rii, a = <^ ia Ci.<j. The coefficient t is 
the hopping energy between adjacent sites, and U is the 
on-site interaction energy. Jaksch et al. (3[ show how a 
gas of atoms in an optical lattice reduces to this model 
at low temparatures, and provided expressions for t and 
U in terms of atomic properties. We consider the case 
U, t > 0. 
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At half-filling (one atom per lattice site) and strong 
repulsion (U 3> t), the Hamiltonian reduces to an anti- 
ferromagnetic Hciscnbcrg model: 

Hafm = J 2^ ^* ' ^J' (•"•) 

with J = At 2 /U and Sj is the spin operator for site i. 
At lower lattice depths, J is renormalized to a slightly 
lower value by the direct interaction between atoms on 
neighboring sites An infinite system described by 
Hafm undergoes a second order phase transition to an 
antifcrromagnetically ordered Neel state at a tempera- 
ture T N ~ J < U. 

Achieving these temperatures is far from trivial. Ex- 
periments on cold fermions 0, [|| have reached suffi- 
ciently low temperatures to observe characteristics of 
the Mott insulator, but so far T > Tjv- Much lower 
temperatures will be reached in the near future, when 
the next generation of cooling protocols are implemented 
0, B B G3 ■ Many of these approaches are interesting not 
only for their utility, but also for the insight they provide 
about fundamental issues of thermodynamics and quan- 
tum statistics. 

The proposals for detecting antiferromagnetic ordering 
in an atomic gas mostly rely on the fact that appropri- 
ately tuned light couples differently to the two spin com- 
ponents [Til . [l3 |. For example, in an experiment where 
one can directly image individual lattice sites [HiHU, one 
can use this selectivity to directly resolve the antiferro- 
magnetic ordering. An alternative proposal involves mea- 
suring the spatial noise correlations in images of the den- 
sity profile after turning off the trapping and lattice po- 
tentials [HI, [T3l • The experimental noise contains in- 
formation about the density-density correlation functions 
of the system, and in principle encodes the antiferromag- 
netism. This method has been used to detect imposed 
lattice order in a bosonic Mott insulator [lj| [TH and a 
fermionic band insulator [2(|. Its disadvantage is sensi- 
tivity to technical noise and that averaging over many 
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experimental shots under similar conditions is needed to 
obtain sufficient statistics. 

Here we propose detecting magnetic order through 
in-situ Bragg diffraction of light (as opposed to Brag g 
diffraction of atoms, another common techniqu e 12 ill) . 
This method is analogous to neutron scattering [22|, [23[ 
and magnetic x-ray scattering [24], [25| , as used in solid 
state physics. Optical Bragg diffraction has been used 
previously in cold atomic gases to confirm the crys- 
talline ordering of nondcgcnerate atoms in optical lattices 
[IH, H3, HH ■ Compared to the other proposed techniques, 
the primary advantages of Bragg diffraction are its sim- 
plicity, speed, and relatively large signal. In this work, 
we report a detailed theoretical analysis of Bragg scatter- 
ing from an array of atoms trapped in an optical lattice 
and demonstrate its usefulness under typical experimen- 
tal conditions. 



II. BASIC THEORY 
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FIG. 1: (Color online) Real (jj — 1, solid red line) and imagi- 
nary (re, dashed blue line) parts of the index of refraction for a 
uniform gas of 6 Li atoms as a function of laser detuning. The 
atoms are taken to be in the two lowest hyperfine states with 
densities n l = n T = l/[2(532 nm) 3 ] = 3.3 x 10 12 cm~ 3 in an 
external magnetic field of 834 G, yielding a Zeeman splitting 
of 2Ao = 2ty x 76 MHz. The laser excites the D2 line with 
linewidth T = 2tt x 5.9 MHz « 2A /13. Thin vertical lines 
denote the two resonances and their midpoint. 



We consider atoms confined in a three dimensional 
simple cubic optical lattice potential, V(x, y, z) = 
Vb [sin 2 (nx/a) + sm 2 (ny/a) + sin 2 (nz/a)], where the lat- 
tice constant a is half of the lattice laser wavelength and 
Vb is the depth of the optical potential, assumed to be 
large enough that the tight-binding approximation holds. 
Although the transverse Gaussian intensity profile of the 
laser beams causes Vq to be position dependent, in the 
relevant part of the cloud it can be taken to be constant. 

Two collisionally closed atomic hyperfine sublcvels (|, 
|) may be treated as a pseudospin-i system, where there 
is no mechanism for spin relaxation. We assume that the 
two states are separated in energy by a splitting 2Ao, 
and have optical resonance transitions suitable for imag- 
ing with line width T - concrete examples using 6 Li will 
be given later. An external magnetic field is applied to 
tune the interactions between atoms via a Feshbach res- 
onance. This field also defines the quantization axis for 
the Zeeman sublevels of the hyperfine states and, hence, 
of the pseudo-spin states. 



A. Bragg scattering cross section in the Born 
approximation 

We initially consider Bragg scattering in the limit of 
low probe intensity and low optical density. In this limit 
the Born approximation may be used. The more general 
case will be considered in Sec. IIV Al and the appendix. 
Our system consists of atoms in an optical lattice de- 
scribed by operators for position ij and occupation num- 
ber fija. In the Born approximation, the total amplitude 
for elastic scattering reduces to [29| 

***** =E E U^e*"*', (2) 



where the difference between incoming and outgoing 
wavevectors of the photons is lq — kj = K, and 

/ct = ~2fc (e ^< A ' ' Gm)(e ™ ' Gkl ' Al) A a TiT/2 
= Jr;( e k/A f • e m )(e r *„ ■ e ki)Ai )e i5 " sin(<5 CT ), (3) 

is the scattering amplitude for individual spin-cr atoms. 
Here, e^A; is the polarization vector of the incoming 
photon with wavevector ki and polarization Xi (gener- 
ally elliptical), and e m is the polarization which couples 
to the resonant optical transition. The phase shifts 5 a 
are related to the linewidth and detuning by tan(£ CT ) = 
— r/2A a . Note that Eq. © is only valid when the optical 
density of the sample is small compared to one. Optically 
dense samples will be considered in Sec. IIV Al To achieve 
spin dependent scattering the probe light frequency must 
be near an atomic resonance. Hence the interaction be- 
tween light and the atoms is neither purely diffractive 
nor purely absorbtive. This is clearly illustrated by con- 
sidering the lattice free case and assuming the atoms are 
uniformly distributed with densities n-\ and n-i . The elec- 
tric susceptibility Xbulk of the bulk sample is then [2^| , 



where to is the angular frequency of the light, c is the 
speed of light in vacuum. The real and imaginary com- 
ponents of the index of refraction 77 and k are then given 
by the usual relationship: 77 + in = y/1 + x- The depen- 
dence of these quantities on the probe laser frequency for 
6 Li is show in Fig. [T] in the case where = n^. 

As seen in the figure, when the probe detuning is 
halfway between the levels, r\ is unity and k has a local 
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minimum. For this detuning the scattering phase shifts 
of the two states obey <5f = —5±. This will be the op- 
crating point of our proposed measurement. There arc 
three advantages to working at this point: (1) we can 
neglect bulk dispersive effects because r\ = 1, (2) absorp- 
tion is small, and (3) we are maximally sensitive to the 
spin order. 

Experiments measure the differential scattering cross- 
section 



da(K) 



= Tr 

9 
Ik- 



A 



<r,<r',j,j' 

xe i ^- 5 -' ) sin(<5 ff )sin((J o .0 



The density matrix pi = pj4<E)|kj, A,) (k,-, A*| describes the 
initial state of the system (where pa is the thermal den- 
sity matrix of the atoms). Tracing over outgoing polar- 
izations gives a factor J2\ f \( e k f ,\ f ' e ™)( e m ' e ki,A«)| 2 = 
i(l + cos 2 6»)(l + cos 2 6>'), where 9 (#') is the angle be- 
tween the magnetization axis and the incoming (outgo- 
ing) wavevector, assuming the incomingelliptical polar- 
ization is chosen to maximize the signal [39j . The thermal 
~ Ti 'n factorizes with (nj a hy a i) = 

zj I \2 1 " ; / 1 where we assume half-filling 
?i; T + riii = 1, cr = ±1, and S zi = i(n lT - n,-;)). The 
factor ( e «(*;-*y)) » e ~ 2W e «(R,-R 3 .0 ; where e -2W = 

is the Debye- Waller factor due to the zero-point 



((l+^i) (3+^'))' 



average \ nj a riji a ie 
1 + 



-l 2 K 2 /2 



motion of the atoms around the lattice sites R, in terms 



of the harmonic oscillator length Z w y^h/muj of the in- 
dividual wells. The differential cross-section becomes 

^ = ii (1+cos20)(1+cos2 ^ 2H/ 



(a(<5 T ,,5 i )C(K)+/3( ( 5 T , ( 5 1 )5(K)), (4) 



where wc introduce the crystal structure factor 
C(K) = Ylij e lK '( Ri-R .») and spin structure factor 
•5( K ) = V, ,«'•<•■"' "i S zi S zj . The coefficients are 

«(*T.*J = ll/t + /il 2 and /?(<Ml) - l/T ~ h\ 2 with 
/ CT = e* 5 " sin ((5,7 ). As described previously, for a laser fre- 
quency where <5j = —5± = 5 <C 1 (this is achieved when 
A 1 = — Aj = Aq), the coefficients simplify: 



sin 4 ((5) r 
sin 2 (2(5) 



5 4 » (r/2A ) 4 , 
4(5 2 w 4(r/2A ) 2 . 



The factor S(K) in Eq. ((4]) represents the well-known 
result that (spin-selective) Bragg scattering measures the 
spin structure factor. In solids one uses neutron scat- 
tering for the same purpose and the coupling between 



the magnetic moment of the neutrons and electron spin 
causes different scattering amplitudes for | , j electrons 
[30lj . In cold atoms the situation is even more favorable 
as this result shows: by choosing detuning with opposite 
phase shifts for the two spin states one can basically turn 
off the signal from the underlying lattice and probe only 
the spin order (i.e. a <C (3). 

The predicted locations of Bragg diffraction peaks 
(i.e. maxima in da(K.)/dQ) are given by the Laue condi- 
tion: the difference in wavevector between the outgoing 
and incoming scattered light (kj — kf = K) is equal to 
a reciprocal lattice vector Q of the crystal [3l| . Rewrit- 
ing this condition in terms of incoming and outgoing unit 
propagation vectors k i; k/, resp. and the wavelength of 
the scattered light A: 



hf — ki = — (I m n) , 



(5) 



where the Miller indices of the scattering plane are de- 
fined for a cubic lattice by (Iran) = Qa/2w. 

Additionally, the magnitude of the reciprocal lattice 
vectors available for Bragg scattering is further restricted 
by the the vector triangle inequality: 



Ik/I-Nl < |Q| < |k,| 
< |Q| < 2k, 



(6) 



where the second line assumes elastic scattering. The 
consequence of the above expression is that for the con- 
figurations typical of optical lattice experiments where 
a ~ A only a handful of choices for K will produce Bragg 
scattering. 



B. Finite simple cubic lattice 

The structure factors S(K) and C(K) in Eq. are 
sums that may be performed exactly for a simple cubic 
lattice with length per side of L sites (N = L 3 total sites). 
By recognizing that the sum in the crystal structure fac- 
tor is that of a geometric sequence in this configuration, 
it evaluates to 



c(k)= n 



sm 2 (K j aL/2) 
sin 2 (A>/2) ' 



(7) 



If K corresponds to a reciprocal lattice vector of the crys- 
tal then the sum is maximized with C(K) = A 2 , demon- 
strating the coherent enhancement of the scattering when 
the Laue condition is satisfied. As K moves away from 
a reciprocal lattice vector, C(K) falls off approximately 
quadratically. 

The spin structure factor S(K.) may be evaluated in 
a similar way. Magnetically ordered phases may be de- 
scribed by a spin ordering vector q such that the average 
staggered spin is s = Se ?q Rj , having z-component 

s z . In the mean field approximation the spin dependent 
factor in 5(K) is (S M S tj ) = (S zl )(S ZJ ) 



s 2 e iq-(Ri-R 3 - 
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In particular, a ferromagnetic (fully polarized) phase has 
q = and a Neel AFM has q = ± ±). With this 

notation 
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s(K)=si n 



sin 2 [(/\j + qj)aL/2] 
siv?[{K 3 + qj )a/2] ' 



(8) 



When K + q corresponds to a reciprocal lattice vector, 
5(K) is maximized with value s 2 N 2 . 

The locations of the maxima in C(K) and S{K) are 
well-known results from solid state physics. As explained 
in standard textbooks [10,[3l|, Bragg scattering for a sim- 
ple cubic lattice occurs only at planes where the Miller 
indices are integers in the absence of magnetic ordering. 
If the crystal contains antifcrromagnctic ordering, half- 
integer Miller indices are also possible, indicating a dou- 
bling of the crystal unit cell. Planes with mixtures of 
integer and half-integer Miller indices have zero scatter- 
ing amplitude. 

In the case of a paramagnetic phase (e.g. when T > 
T/v), the spins S{ Z appearing in the sum S(K) do not 
show long range order. We can then estimate the dif- 
fuse background by assuming that all the spins point in 

random directions: 5(K) = (£V 5«e iI£ - Ri > 2 « A/4, in- 
dependent of K, i.e. this term describes incoherent scat- 
tering. For a real experiment at finite temperature, a 
portion of the sample is likely to be in the paramagnetic 
phase. The incoherent scattering from these disordered 
spins will yield a diffuse background. 



III. 



EXAMPLE: 6 Li IN A RED-DETUNED 
LATTICE 



To demonstrate this technique for a typical experi- 
mental configuration, we simulate a crystal of 6 Li atoms 
cvaporativcly cooled to Fermi degeneracy in a far off- 
resonance optical trap, then loaded into a 3d simple- 
cubic lattice created from three orthogonal far red de- 
tuned laser plane- wave standing waves at 1064 nm. This 
experiment is currently underway at Rice University. 

The 6 Li atoms are prepared in an incoherent balanced 
mixture of the F = i, nip = ±-| Zeeman sublevels of the 
2Si atomic ground state (denoted J, and |, respectively). 
These two sublevels play a role analogous to the spin- 
^ electrons in crystalline solid. An externally applied 
uniform magnetic field controls interatomic interactions 
via a Feshbach resonance (H)]. We choose the direction 
of the magnetic field to be along one of the lattice axes 
for experimental convenience. This field also determines 
the quantization axis of the atomic states. 

The probe light is tuned near the 2Si — > 2P| reso- 
nance of 6 Li with wavelength A = 671 nm. Because A is 
similar in magnitude to the lattice constant a, only a few 
lattice planes will produce diffraction. Applying Eq. ©: 
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FIG. 2: (Color online) Geometry for scattering from (0 1) 
plane (left) and (555) plane (right) for a cubic lattice with 
lattice spacing a. The scattering planes are shaded. The thin 
solid line is normal to the scattering plane. The color of the 
lattice sites indicates spin; this figure shows AFM ordering. 
The arrows represent the directions of the incoming and out- 
going scattered light. The labeled angles are between the out- 
going light and the normal, which is the same value as the an- 
gle between the incoming light and the normal. The values of 
the angles are given by the Bragg condition: cos _1 (|Q|/2|k|), 
assuming lattice spacing of a — 532 nm and probe wavelength 
of 671 nm. 



The only valid nonmagnetic scattering planes for our 
system are (±1,0,0), (±1,±1,0) and permutations of 
these indices. (The (0 0) case represents unscattered 
light). The allowed magnetic scattering planes are 
(± i , ± i , ± i ) . The basic geometry for two of the scat- 
tering planes is depicted in Fig. [21 with the values of 
the scattering angles shown. In an experiment one could 
probe a nonmagnetic peak such as (0 1) to confirm the 
cubic ordering of the sample and a magnetic peak such 
as (| I |) to detect the AFM ordering. The ratio of 
these two peak intensities is a measure of the staggered 
magnetization, proportional to s 2 and independent of N. 

The Zeeman splitting between states | and J. is 2Aq = 
2tt x 76 MHz m 13T when the external magnetic field is 
tuned near the Feshbach resonance at 834 G. The bulk 
optical properties of the sample are those shown earlier in 
Fig. [U for one atom per site and equal spin populations. 

To estimate the optical signal strength produced by the 
scattering we consider a detector (CCD array or photo- 
diode) placed r = 0.5 m away from the atoms along the 
direction of the outgoing wavevector where we expect a 
peak. We then numerically evaluate Eq. ([4]), assuming 
a perfectly ordered sample of N = 41 3 atoms and a lat- 
tice depth of Vq = 10 /iK (about 7 lattice photon recoil 
energies). The value of the spin structure factor varies 
depending on the phase we wish to calculate. For a fully 
polarized phase the spin ordering vector is chosen to be 
q = 0. For the AFM phase, we use q = 2zt(I 1 1). 
Finally, for the paramagnetic state we assume the spin 
structure factor is independent of direction with a value 
of S(K) = N/A, as mentioned previously in Sec. Ill Bl 
The maximum light intensity for the Bragg peak is given 

by /out (k' ,?") = % dcr j^ - , where |£ denotes the differen- 
tial cross-section in direction k'. We take the probe laser 
to be a plane wave with intensity I la = 0.5 mW cm' 2 
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TABLE I: Predicted experimental values of maximum scat- 
tering cross-section, peak intensity, and scattered power for 
various magnetic states (paramagnet [PM], antiferromagnet 
[AFM] and polarized [Pol.]) and scattering planes using typ- 
ical parameters, calculated as described in the text 
PM AFM Pol. 



m K) 

I (W/cm 2 ) 
P (W) 

& ( cm2 ) 
I (W/cm 2 ) 

P (W) 



(0 1) scattering plane 
1.4 x 10" 5 1.4 x 10" 5 2.3 x 1CT 3 
2.7 x 1(T 12 2.7 x KT 12 4.5 x 10" 10 
5.7 x 10" 12 5.7 x 10" 12 9.4 x 10" 10 

(|||) scattering plane 
4.7 x 1(T 8 3.2 x 1(T 3 7.3 x 10 
9.4 x 1(T 15 6.5 x KT 10 1.5 x 10" 19 
.0 x 1(T 14 1.4 x 1(T 9 9.3 x 1(T 19 
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and beam waist much larger than the sample size. The 
total scattered power is estimated by numerically inte- 
grating over the area of the detector from the center of 
the peak radially to the first minimum. The results of the 
calculations for the various phases and the two scattering 
planes of interested are given in TablcUl For the maximal 
scattering configurations wc find approximately 1 nW of 
power in the scattered beam spread over about lOmrad 
of angular divergence. Using commerically available pho- 
todiodes (with an appropriate collecting lens) and inte- 
gration times of 10 /xs should be sufficient to detect this 
power. For higher sensitivity, inexpensive CCD cameras 
may be used [40(. 

By varying the direction of the incoming beam, we find 
that samples of this size are sensitive to the direction of 
the incoming beam on the ~10mrad scale. Given the 
size of our simulated crystal, this is consistent with the 
expected diffraction limit. 



A. Finite temperature 

The previous analysis assumed that the sample was 
completely ordered. At finite temperature, especially for 
T near Tn, the AFM ordering will be incomplete. To the 
estimate effect of finite temperature on the Bragg scat- 
tering intensity, we evaluate Eqs. ([J|, |[7J), and ||SJ) for a 
system with L = 41, assuming the staggered magnetiza- 
tion as calculated by a quantum Monte-Carlo calculation 
for a smaller system with L = 6. While providing an es- 
timate of the role of thermal fluctuations, our procedure 
of scaling up the 216 site result to 68921 sites will over- 
estimate the finite-size smearing of the phase boundary. 
We used the ALPS (Algorithms and Libraries for Physi- 
cal Simulations) package to calculate this magnetization 
[33l ] . As shown in Fig. [3] this extrapolation predicts that 
the (522) B ra §S peak begins to appear as the tempera- 
ture goes below Tjy = 0.95 J, which is consistent with the 
bulk behavior of the system [34j ■ By measuring the ratio 
of the two Bragg peak intensities we gain a measurement 
of s z independent of shot-to-shot fluctuations in N. 



FIG. 3: (Color online) Staggered spin per atom of the system 
(red squares) calculated using a quantum Monte Carlo simu- 
lation of 6 3 atoms, and maximum scattering cross-sections of 
the magnetic and nonmagnetic Bragg peaks (blue circles and 
green triangles, resp.) as a function of temperature. The 
cross-sections are calculated using Eq. Q with the above 
value of staggered magnetization for a system of 41 3 atoms. 



IV. FURTHER CONSIDERATIONS 

A. Multiple Scattering 

For optically dense samples multiple scattering can- 
not be neglected when calculating the differential cross- 
section. We calculate the effect of multiple scattering by 
numerically solving the Lippmann-Schwinger equation, 
which is derived in the appendix. The key result is the 
scattering amplitude (also Eq. (|A.3I) ) 



Ft 



3 hT , 
~2fcT ^ ek/ ' A/ 



(9) 



where Aj are numerical coefficients found by solving 
Eq. CO]). 

Using these results, we qualitatively study the con- 
sequences of finite optical density (OD), in particular 
the interplay with the position uncertainty of the atoms. 
(The Debye- Waller factor in Eq. (H|) is an approxima- 
tion of this effect [Hj.) In this section we consider a 
spherical sample and assume that all lattice sites within 
a radius R are occupied. Position uncertainty caused by 
zero-point motion of the atoms is modeled as random 
displacements of the (point-like) atoms by a Gaussian 
distribution with standard deviation a = Ax x a. For 
the ground state of a sinusoidal well in the tight-binding 
regime Aa; = 1/ttVq^ 4 , where Vq is the lattice depth in 
units of the lattice recoil energy. For a typical experi- 
mental value of Vo = 7, Aa; = 0.2. The lattice spacing 
is a = 532 nm and the wavelength of the probe beam is 
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FIG. 4: (Color online) (a-c) : Differential cross-section of the 
(0 1) Bragg peak for a fully polarized spherical sample with 
N = 925 atoms (OD = 37), resonant probe light, and vary- 
ing mean displacements (disorder) Ax = a/a = 0, 0.1, 0.2 
(corresponding to (a), (b), (c)) calculated from Eqs. ([9| and 
(|A. 1|) (blue triangles), as a function of outgoing direction 59 
relative to the peak center (directions shown in Fig. [2]). For 
comparison we also show the Born approximation result (red 
circles). The lines are guides to the eye. (d): Intensity of the 
(0 1) Bragg peak as a function of detuning A in the fully 
polarized case Ax — 0, 0.1, 0.2 (solid red, dashed green, dot- 
ted blue lines). The resonance is at A = 0. (e-f): Intensity 
of the (5 | |) Bragg peak as a function of A for a Neel AFM 
state with Zeeman splitting between spin states of 2Ao = 2r 
for varying Ax = 0, 0.1, 0.2 (solid red, dashed green, dotted 
blue lines) and particle numbers (e): N = 925 (OD = 7.4 at 
A = 0) and (f): N = 257 (OD = 4.8 at A = 0). The data 
in (d-f) have been symmetrized about A = to remove the 
effect of residual polarization caused by the small sample size. 
The vertical dotted lines in (e, f) indicate the atomic optical 
resonance. 



A = 671 nm. 

Fully polarized: We first consider a model for a fully 
polarized state (all atoms in the same pseudospin state) 
with R = 6a (N = 925 atoms). The angular distribution 
of the differential cross-section is shown in Fig. 0Ja-c) 
for various Ax using both multiple-scattering and Born 
approximation methods. The width of the peak is con- 
sistent with the diffraction limit for a sample of this size. 
As shown in Fig.QJd), we find that for high optical den- 
sity, the intensity of the (0 1) Bragg peak develops a 



local minimum at resonance (zero detuning). On res- 
onance the optical density is highest, and even though 
the optical scattering cross-section for individual atoms 
is maximal, the Bragg signal has a local minimum be- 
cause of absorption. This effect is similar to that seen in 
the experiment by Birkl ct al. (26|, except in that case 
the disorder was caused by temperature. 

AFM: The antiferromagnetic cloud is modeled by as- 
signing lattice sites different pseudospin states according 
to Neel order. We define A as the laser frequency detun- 
ing relative to the center of the two atomic resonances. 
We assume a small Zeeman splitting Aq = T in the sim- 
ulation so that the small sample of 925 atoms is optically 
dense when the probe laser is tuned directly between the 
resonances (OD = 7.4). This splitting is much smaller 
than the experimental parameters for 6 Li near the Fesh- 
bach resonance, where Ao ~ 6.5T. The optical density of 
the simulation is comparable to what one expects with 
~4x 10 7 6 Li atoms [4l| . about 40 times the sample size 
expected in the experiment. 

The results are shown for different radii R (i.e. different 
particle numbers) and disorder Ax in Fig.[4^e, f). We find 
that the maximum intensity of the {\ \ \) Bragg peak 
scattered from an optically dense sample (e.g. Fig.0Je)) 
is not found on one of the resonances (as one finds in 
the Born approximation or samples with low optical den- 
sity), but rather around the midpoint in between the | , [ 
states (A ss 0). Disorder increases incoherent scattering, 
reducing the amount of light coherently scattered into 
the Bragg peaks - the disorder effectively enhances the 
optical density of the sample in a nonlinear way. 

Under typical experimental conditions (finite temper- 
ature in a harmonic trap) the outer shell of the atomic 
cloud will be disordered and we expect the AFM order to 
develop only at the center. The disordered atoms could 
then shield the AFM from the probe light if the optical 
density in the shell becomes too high. The above calcu- 
lation demonstrates that in this case it is favorable to de- 
tune away from the t, I states to increase the penetration 
depth of the light and maximize the signal at the (h h | ) 
Bragg peak. The physical picture is that for sufficiently 
large detuning (but between the two resonances) the dis- 
ordered portion of the sample only contributes weak in- 
coherent scattering, but the AFM portion still exhibits 
detectable coherent Bragg scattering from the magnetic 
planes. 



B. Symmetry effects 

The ordered state of a classical antiferromagnet is char- 
acterized by having a nonzero staggered magnetization s. 
In the absence of a symmetry breaking field, s will point 
in a random direction on the Bloch sphere (see Fig. [5]), 
and sequential experiments will find this order pointing 
in a new direction. Quantum mechanically the physics is 
similar, with the added feature that the staggered mag- 
netization may be in a coherent superposition pointing in 
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FIG. 5: (Color online) Cartoon of the Bloch sphere repre- 
sentations of the various AFM configurations. Left: isotropic 
AFM state, where the spins (large red arrows) may lie any- 
where on the Bloch sphere. The z-component of the stag- 
gered spin vector is shown by the small blue arrow. Center: 
canted AFM state, where the spins are restricted to be near 
the xy-p\ane but lift off of it by some small canting angle (ex- 
aggerated here for clarity). In this case the z-projection of the 
staggered spin vector is zero. Right: canted AFM state after 
applying a tt/2- pulse. Now the staggered spin has a nonzero 
z-component (small blue arrow), which may be detected by 
Bragg scattering. 



all directions. From a practical standpoint, this distinc- 
tion is minimal: a measurement will yield a definite (and 
random) direction for s. This (classical or quantum) un- 
certainty means that the the scattered intensity J of a 
magnetic Bragg peak will fluctuate from experiment to 
experiment. The scattered light only interacts with the 
t and I eigenstates, and hence is only sensitive to the 
z-component of the order s z . Assuming an isotropically 
distributed s, the probability distribution of the inten- 
sity will be ff (I) = \//max/4/, with a mean value of 
iiso = /max/3, where 7 max is the intensity when the spins 
are parallel to z. This purely geometric effect necessitates 
repeating the experiment multiple times. 

Recent numerical calculations point out that a suffi- 
cient polarization of the sample (|fif — njj > 0) favors 
a staggered magnetization lying in the xy-plane of the 
Bloch sphere [III [36J. This response to a spin imbal- 
ance is analogous to the response of an electronic an- 
tifcrromagnct to a magnetic field and can be explained 
in terms of a n anisotropic spin susceptibility. As il- 
lustrated in Fig. [5j it is relatively easy to add a finite 
polarization to xy-antiferromagnetic order: one simply 
cants each spin slightly towards the z-axis of the Bloch 
sphere. On the other hand, it is much more energeti- 
cally costly to polarize a z-antiferromagnet: one must 
completely flip individual spins, adding 2 J energy to the 
system for each spin that is flipped. Thus, the lowest 
energy configuration of a polarized AFM has the Neel 
order in the xy-plane. One expects that even at finite 
temperature the thermodynamic ensemble will be domi- 
nated by configurations with s lying in the xy-planc. A 
simple model is to take the spins of the two sublattices 
to be Si, 2 = \ [±(cos</>cos£)i ± (sin</>cos£)y + (sin£)z], 
where 4> is uniformly distributed over 2tt radians of a cir- 
cle and £ <C 1 is a small fixed canting angle relative to 
the xy-plane, which is set by the polarization. 

In the canted configuration, the z-component of the 




FIG. 6: (Color online) Shot-to-shot probability distributions 
for the intensity of the Bragg scattering signal for two cases. 
The solid red line shows the case when the staggered mag- 
netization may point in any direction. The dashed blue line 
shows the case for a canted AFM, after applying a 7r/2-pulse 
to rotate the spins from the xy-plane to the rrz-plane. The 
vertical lines indicate the mean values of the respective distri- 
butions. The vertical scale is normalized such that the areas 
under both curves are one. 



staggered magnetization is zero (center frame of Fig. [5j , 
making the AFM invisible to Bragg scattering or opti- 
cal imaging. However, the spins may be reoriented to 
partially project along the z-axis by an RF 7r/2-pulse at 
the Zeeman splitting frequency, giving a non-zero Bragg 
signal (right panel of Fig. [5]). The resulting shot-to-shot 
intensity probability distribution of the Bragg scattered 
light is then ff™\I) = - J// max ) with 

a mean of / C ant = -fmax/2. The intensity distributions for 
Bragg scattering of the isotropic AFM and canted AFM 
(after applying the 7r/2-pulse) are shown in Fig. [5] Note 
that the 7r/2-pulse does not change the probability distri- 
bution of the isotropic case because of symmetry. Mea- 
suring these probability distributions provides a method 
to experimentally distinguish between the isotropic and 
canted AFM states, particularly because of the different 
behavior of the two distributions when / ~ /max- 



V. CONCLUSION 

These calculations confirm that optical Bragg diffrac- 
tion is a viable method to observe antiferromagnetic or- 
dering of atoms in an optical lattice. The experimental 
measurement consists of using a CCD or photodiode to 
measure the intensities of the peaks. The presence of 
the (0 1) Bragg peak indicates cubic ordering, while 
the (555) Bragg peak indicates antiferromagnetic or- 
dering. In the actual experiment, only a portion of the 
atoms near the center of the crystal will be in the anti- 
ferromagnetic state. The ratio of the (| | 5) intensity 
to the (0 1) intensity provides a measure of the degree 
of ordering. Additionally, the shot-to-shot intensity dis- 
tribution provides a method to determine whether the 
antiferromagnetic state is isotropic or canted. 
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APPENDIX: FULL SCATTERING FORMALISM dimensionless function Qj X = /3 jt - i-yji has [3i 



Here we derive the full transition matrix (T-matrix), 
including multiple scattering events. We consider scat- 
tering of a photon off of TV 2-level atoms located at posi- 
tions Yj described by the matter-light Hamiltonian 

H = huj 0j S z j + ^2 ^k.AO^Ok.A ~y^dj ■ E(rj), 

3 k.A j 



Ho 



where Tiu)o t j is the energy splitting between the two levels, 
S z .j = (\e)j {e\j - \g)j (gU /2 is the z-componcnt spin 
operator for the j-th atom in terms of the corresponding 
ground state |<?) • and excited state |e) •, w k .A is the fre- 
quency of a photon with wavevector k and polarization 
A, \ and at. a are creation and annihilation operators 
of the same photon, d = er is the dipolc operator (e is 
the electron charge) and E(r) is the quantized electric 
field at the location of the atom at r. In the rotating 
wave approximation (RWA) the matter-light interaction 
becomes [I?} 

V = - i 3k [e ikrj (e k ,A ■ e; n )S+a KX 

k,A,j 



where De*^ 
S+ = \e) (g\, Sj 



(e\jdj 



Iff) 

\g)( 



e = z, e± 
and g k 



T(x ± iy)/ V2, 
D v / 2nhuj k /L 3 (L 3 



denotes the volume of the sample). In the low inten- 
sity limit (much less than the saturation intensity of 
the atomic transition) we may consider scattering of a 
single photon. Then within the RWA the Lippmann- 
Schwinger equation for the outgoing photon = 

\ki,Xi) + (E - H + te^V |\p(+>) (where |k;,A;) is an 
incoming plane wave photon with wavevector and po- 
larization Xi , e is a positive infinitesimal) closes with the 
ansatz 



*(+>)= £>k,A|k, A)+J^ViS 



"iji 
where r 



cos(krji) 



kr 



ji 



sm(krji) cos(krji) 



(kr 3l y 



{kruf 



sm(krji) 
kr 3 



(krj 

ri\ and k = |k,- 



cos(krji) sm(krji) 



■jl = \rj — aim k = = \kf\ [42fl . and p and 
q depend on whether one considers a m = or m = ±1 
transition. One has p = 1 — (z ■ rji) 2 , q = 1 — 3(z • iji) 2 
for m = and p = |(1 + (I ■ r^/) 2 ), g = |(3(i • r jZ ) 2 - 1) 
for m = ±1 I38H. The transition matrix is 



/ k . . A . : k A = (kf,\f\V\^ +) ) 
2 



9k ( e k/,A f • <,) (e m • e kijAj 
' k r ,1,. 

i 



(A.2) 



Eq. (|A.1[) includes the effects of multiple scattering to all 
orders, provided one can solve for Aj. In the limit of low 
optical density of the sample, multiple scattering can be 
neglected and one obtains the Born approximation result 



Aj = (AEj 1 e lW > r i or 



I'k. x. k ..\ =g k {ek f ,\ f ■ e* n ) (e r „ • e k . 



J2e^-^(AEj+i'^- 



We finally note that the scattering amplitude Fk.\. jk^Aj 
is related to the T-Matrix through 



k f ,A/;ki,Ai 



_3_nr . , . 

2k 2 ^ 6k/ ' A/ ' Gm > 



k.A 



,k ' '• ,1,. (A.3) 

3 

In the Born approximation this reduces to Eq. 



and one obtains a linear equation for effective amplitudes 
Aj (related to tpj by = -ig k (el j X . ■ e m ) Aj) 

(AEj + i^fj A, = + ^ £ 9**' (A' 1 ) 



where AEj is the (renormalized) detuning at site j and 



T = Auj'^D 2 /(3hc 3 ) is the linewidth of the transition. The 
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